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ABSTRACT 

Previous results have D>2,3], been concerned 
with establishing (1) whether operating problems 
could exist when WECS generation Is significant and 
(2) the proper modification of unit commitment, regu- 
latlon, and economic dispatch required to provide 
sufficient system security and alleviate the oper- 
ating problems caused by WECS generation changes. 
This" paper discusses methods of reducing the WECS 
generaTfon change through selection of the wind tur- 
bine' ri^daT for each site, selection of an appropriate 
-Siting configuration, and wind array controls. An 
analysis of wind generation change from an echelon 
and a farm for passage of a thunderstorm is presented 
to estjblish the factors concerning the wind turbine 
model and siting configuration that contribute to 
these variations. Detailed simulation results indi- 
cate more precisely how these factors can be 
exploited to minimize the WECS generation changes 
Reduction of the wind generation change 
ovBT ten minutes Is shown to reduce the Increase in 
spinning reserve, unloadable generation and load fol- 
lowing requirements on unit cormitment when signifi- 
cant WECS generation is present and the farm 
penetration constraint Is satisfied. Controls on the 
blade pitch angle of all wind turbines in an array or 
a battery control are shown to reduce both the wind 
generation change out of an array and the effective 
farm penetration in antictpation of a storm so that 
the farm penetration constraint may be satisfied. 

. r::;::. ; i. introduction 

The research reported in this paper is an exten- 
sion c-f previou$~work [1,2,3]. The objectives of the 
earl 1 er \^rk "was to determi ne : 

{!) if operating problems could exist on auto- 

_ ; matlc generation control (regulation and 

economic dispatch), frequency regulation, 
and unit conmltment when wind generation 
capacity Is significant; 

(2) the penetration limits on wind generation 
capacity that would alleviate these oper- 
ating problems. 

The results indicated that there were two spe- 
cific operating problems which could both be elimi- 
nated with proper penetration constraints: 

(1) The automatic generation control will satu- 
rate for long periods when the total change 
in wind generation for passage of a thun- 
derstorm front and simultaneous load change 
in a ten minute interval will require non- 
wind generation change that exceeds load 
following capability in a ten minute inter- 
val. This problem violates NAPSIC perfor- 
mance standards but can be eliminated by 
imposing a farm penetration constraint on 
the capacity of all wind turbine generators 
that can be affected by a single thunder- 
storm front. 

(2) A cycling problem caused by simultaneous 
load and generation change that Induce fre- 
quency deviations that exceed governor 
deadband. This continual cycling of steam 
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turbine units is objectionable to generator 
operators and can cause increased main- 
talnance costs, forced outage rates and 
ultimately reduce unit life. The cycling 
of nuclear units is of concern for safety 
reasons in addition to those mentioned 
above. The cycling problem can occur due 
to a storm front sweeping through a wind 
generator array causing large power varia- 
tions on successive echelons. A echelon 
penetration constraint on the capacity of 
all WTGs that can experience simultaneous 
change In generation level will eliminate 
this cycling problem. 

A subsequent study [3] was devoted to a detailed 
discussion of the modification of unit commitment, 
regulation, and economic dispatch when WECS genera- 
tion Is significant, A modified farm penetration 
constraint is determined that limits WECS generation 
to be less than the maximum first contingency loss of 
resource or commitment. A violation of this farm 
penetration constraint is shown to necessarily cause 
an Increase in the maximum first contingency loss of 
resource or commitment to the level of the farm 
capacity and thus an increase In load following, 
spinning reserve, and unloadable generation require- 
ments on unit commitment, A discussion of the 
methodology, costs and benefits of changing unit 
commitment, when WECS generation Is significant and 
the farm penetration constraint is or is not vio- 
lated, is included. A discussion of the methods for 
modifying unit commitment is also included. Detailed 
simulation results that document the reduction of the 
effects of significant WECS generation change through 
the modification of the unit commitment regulation, 
and economic dispatch is also presented. 

A modified echelon penetration constraint is 
proposed which limits instantaneous change and rate 
of change of wind array generation that must be 
handled by frequency regulation and regulation con- 
trols. This constraint is imposed to limit cycling 
of units. 

This paper presents both analysis and simulation 
results that show how to decrease the WECS generation 
change over a ten minute interval through selection 
of the wind turbine generator model at each site, the 
siting configuration, and controls on the power 
variation out of the array. These direct controls of 
power out of the array are shown to also permit 
reduction of the effective farm penetration below the 
farm penetration constraint level and thus make the 
increase in spinning reserve, unloadable generation, 
and load following requirement depend on the probable 
change in WECS generation over a ten minute Interval 
rather than on the farm capacity, which would be the 
maximum first contingency loss of resource or commit- 
ment for a particular utility if the farm penetration 
constraint were violated. 


2. ANALYSIS AND SIMULATION OF WIND GENERATION 
CHANGE FROM AN ECHELON AND A FARM 

The purpose of this section is to: 
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(1) briefly describe the model and simulation of 
wind power change from an array of wind 
turbine generators for passage of a thun- 
derstorm front; 

(2) analyze the power change and rate of change 
from an echelon and a farm in terms of the 
factors that determine these changes; 

(3) review the methodology for determining 
spinning reserve, unloadable generation, 
and load following requirements on unit 
commmitment and briefly discuss how the 
selection of wind turbine models, siting 
configuration and wind array controls can 
influence these requirements. 

A model of a single MOD-1 WIG and an array of 
wind turbine generators is developed. The MOD-1 WIG 
model, given in [l], is a static nonlinear model that 
relates generation to wind speed if the wind speed 
does not exceed the cut out velocity for a sustained 
period jWhich causes shutdown to avoid damage. The 
dynamics of the shutdown startup sequence is also 
modeled since a thunderstorm can cause such a shut- 
down. A similar model of a MOD-2 wind turbine is 
discussed in section 3 of this paper in order to 
compare the changes from an identical siting configu- 
ration of MOD-1 and MOD-2 wind turbines experiencing 
the identical wind speed profile. 

The farm model, which is cormnon for every wind 
turbine model, assumes the motion of thunderstorm 
front is normal to each echelon. The power out of the 
first echelon is then Just the generation out of a 
single generator in this echelon multiplied by the 
number of generators in this echelon. The generation 

out of the echelon is the generation profile of 
this WTG in the first echelon (1) delayed by an 
interval (d./V ) proportional to the distance between 
J ®th 

the first and j echelon and inversely proportional 
to the speed of the thunderstorm front, and (2) 

multiplied by the number of WTGs in the echelon. 
The generation out of all echelons is simply sunmed 
to obtain the generation time profile for the farm 
for passage of a thunderstorm front. 

The worst case change and rate of change from a 
coastal farm will now be determined. The results are 
derived based initially on simulation of a worst case 
coastal farm experiencing a worst case thunderstorm 
wind speed versus time profile. 

The worst case MOD-1 siting configuration for a 
coastal farm D]> which is a farm located on the coast 
of a body of water, shows a 0.5 ml. spacing between 
50 generators in each echelon and a 2 mi. separation 
between two echelons. The echelons are assumed sited 
normal to the motion of the thunderstorm front. 

The wind speed profile of a worst case thunder- 
storm gust front on the WTGs in the first echelon is 
shown in Figure lA. The initial wind speed increase 
is due to the leading edge outflow and the second 
segment of high wind speed is due to the trailing edge 
inflow. The wind speed Increases from 13 km/hr. at t 
= 0 to 26, km/hr. at t = 50 causing the power on each 
WTG in the first echelon to increase from zero to 
capacity (1,5 in that interval. The power re- 
mains at capacity for speeds in excess of 26 km/hr, 
due to blade pitch controls. The thunderstorm front 
was chosen with a 13 km/hr. initial wind speed to 
cause maximum power variations out of any WTG. 


The power variation out of the coastal wind farm 
of MOD-1 WTGs, shown in Figure IB, shows two 75 
megawatt ramps each 50 seconds long which are the 
Increases in generation due to the leading edge out- 
flow passing over the two echelons. The two pro- 
nounced power decreases are due to the shutdown of 
WTGs on both echelons caused by sustained wind speed 
exceeding cutout velocity. The time interval between 
the successive Increases or decreases on the two 
echelons is 240 seconds. 
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Figure 1 Simulation of Coastal Farm for a 
Thunderstorm Front 

The logic for initiation of shutdown of a WTG 
requires the output of a one minute smoothed wind 
speed W^y(k), to exceed 64 km/hr. thus, a shutdown 

only occurs for the trailing edge Inflow because the 
excessive wind speed for the leading edge outflow Is 
not sustained long enough to trigger a shutdown. 

The power out of any MOD-1 WTG does not decrease 
after the first 50 seconds of the leading edge out- 
flow passes over (until the shutdown) due to the 
blade pitch control that maintains constant maximum 
generation over a wide range of wind velocities (26 
km/hr, - 64 km/hr.). 

An analysis of power variations out of an eche- 
lon and a farm Is now performed to determine the 
factors that Influence power increases and decreases 
for passage of a thunderstorm front. The analysis 
assumed all WTGs in a farm are similar, and arranged 
in straight parallel lines normal to the motion of 
the front. The analysis is not restricted to any 
particular WTG model if the parameter D is inter- 
preted as the distance between the leading edge of 
the thunderstorm and the point where the wind speed 
reaches Vj^, the wind speed at which that WTG model 

achieves rated generation. This maximum change and 
rate of change in generation Is derived based on the 
additional assumption that the wind speed is below 
V^., the cut in velocity for the WTG before the front 

arrives. 

The time in seconds for a particular WTG to 
change its generation from zero to capacity for a 

thunderstorm front is 

T^ = 3600 D/V^ sec (1) 

where is the velocity of the front and D is the 

distance from the very leading edge of the front to 
the point internal to the front at which wind speed 
first reaches Vj^, the wind speed level just suffi- 
cient for maximum generation on that WTG model (C^). 

Thus a thunderstorm front with a minimum value of T„ 

M 

due to a minimum value of D and a maximum value of V , 

o’ 
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would require a higher response rate for the governor 
frequency regulation and AGC regulation co"ntrd1s to 
handle this change in wind power generation without 
excessive, or sustained chsage in frequency or area 
control error. 



The time Interval between initiation of gen- 
eration changes on two adjacent echelons is 


T- 3600 d 

T, -S-- 


sec 


( 2 ) 


yhere d is the distance between echelons in miles. 
The distance d must be greater than D for the response 
of two adjacent echelons due to passage of the 
leading edge outflow not to overlap. The shorter T^ 

and”'d“ the higher the response rate capability of the 
power" system required to handle this generation 
change without excessive or sustained frequency or 
ar^a control error changes. 


The maximum change of generation from an echelon 
of Ng WTGs with rated' capacity Cy is 


and the"ir^ 3 xjmum average rate of change during period 
required for passage of the leading edge outflow 

is 




(3) 


The maximum average rate of change of power from an 
^^^Selon during the period required for shutdown 

during the passage of the trailing edge inflow of a 
thunderstorm front is 


"APf N 


e 

ht 


■e * 


(4) 


The maximum power change out of a farm composed 
of WTGs for passage of either the leading edge 

outflow or trailing edge inflow Is 

(5) 


aP^ - Nj, X Cy 



The maxT mum" average’ rate of change from the farm 
during tfie“parTodlM^/Ng - 1) Tg + required for the 

fncrease in generation during the passage of the 
leading odge outlfow is 


AP 


1 ^ 
t 


Nf X C„ 




(6) 


e m 


The maximum average rate of change from the farm 
duri ng t^e period (N^/N^ - 1) T^ + T^ required for 

" ^utd owy of ^air^ during passage of the 

trailing edge'Tnflow of the thunderstorm front is 



AP 


At 


Nf X Cy 


fi-) 


(7) 




The formulas for WECS generation change out of 
an echelon and farm for passage of the leading edge 


outflow are given by equations (3) and (6) where T^^ 
and Tg satisfy^ equations (1) and (2) respectively. 

Results obtained In section 4, where generators In 
echelons are randomly sited 1n a strip D miles long 
rather than In straight lines normal to the Ration of 
the front. Indicated that the rate of change of WECS 
generation in an echelon for passage of a thunder- 
storm trailing edge inflow to be identical to that 
for the passage of the thunderstorm's leading edge 
outflow. Thus, the formulas for rate of change of 
power from an echelon (3) or a farm (6) are appropri- 
ate for passage of both the leading edge outflow or 
trailing edge inflow if the siting configuration is 
not in straight lines normal to the motion of the 
front which will generally be true. 

These formulas will not be used to derive de- 
tailed expressions that indicate more precisely the 

dP dP. 

factors that contribute to and so that the 

model WIG and siting configuration can be selected to 
keep these WECS power generation rates below that of 
the power system average response rate capability. 


The power rate of change out' of an echelon is 
quite different depending^ on whether the density of 
wind turbines in an array~Ts of Ts not uniform; l.e.. 


whether the density of wind turbine In an echelqn, 
defined 


^ D X 0_ 


(8) 


is or is not Identical to the density in the farm 


Nf ^ /Nf \ 

" (TT + D) DV 

0 0 ' e * 

where I is the time for the thunderstorm to nrove 
through the farm. The formula for power rate of 
change from an echelon (3) becomes upon substitution 
of (?) and (9) 



Pf Vo Cw 
35DtJ — 


( 10 ) 


if == Pg as it is more nearly in a low density 
midwestern farm siting configuration. Note that 
and Dq are the velocity of the thunderstorm front and 
the width of the front respectively, C^ is the capac- 
ity in MW for the wind turbine model, and is the 

— T 2 

uniform density of wind turbine in the farm in #/m1. 
The rate of change of power out of an echelon is 



p D V 
0 0 

“TdOO 





]) 


(n) 


if pg > Pf, as in the coastal farm siting configura- 

dPe 

tion. Note that jn_ .this .case depend on 




The slope of the power versus wind 
speed curve for a wind turbine 
model for the range of velocities 
Vr - V^j where power can change. 

The slope of the wind speed profile 
of thunderstorm front for the dis- 
tance D into the front. 
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D - The distance Into a particular 

thunderstorm front where power will 
change on the WTG model considered 
or alternately the length of an 
echelon in the direction of motion 
of the thunderstorm front. 

The formula for power out of a farm (6) has the 

form 


dP, 

■at 


^f °o ^0 Cw 
3500 — 


(12) 


(l-k+i - L. )T - predicted change in thermal load In 
ten minutes 

T = .1667 hour = ten minutes 

Q^k Maximum probable drop in WECS genera- 

tion output in 10 minutes 

Q^k Maximum probable rise in WECS genera- 

tion output in 10 minutes 

Q^k Maximum probable rise in system load 

in 10 minutes 


upon subsititution (1,2, 9‘) into (6). Note when the 
density of an echelon is the same as that of the farm, 
the rate of change out of an echelon is the same as 
that out of the farm. 


It Is clear that the wind generation rates of 
change from an echelon (10,11) and farm (12) depend 
on the width (0^), speed (V^^), and the slope 

- \i 

Q Of the leading edge outflow for the thun- 
derstorm front. The remaining parameters depend 

either on the WTG model (C^, ^ — , D) or on the 

siting configuration (p^, , d). The effects of 

different WTG models and siting configuration charac- 
teristics on the change and rates of change from an 
echelon and a farm for different wind speed profiles 
will be demonstrated in the next two sections of the 
paper. The purpose of studying the effects of WTG 
models and siting configuration patterns is to ana- 
lyze how such factors can be used to minimize WECS 
generation change and thus minimize any Increase in 
load following, spinning reserve and unloadable gen- 
eration requirements provided through unit commitment 
modification when WECS generation is present. These 
increased requirements add to fuel, operating, and 
maintainance costs but are required to maintain sys- 
tem security as discussed in [3]. Changes in regula- 
tion and economic dispatch (load following) controls 
must also be implemented [3] to take advantage of the 
increased response and response rate capability pro- 
vided by the increased spinning reserve, unloadable 
generation, and load following capability, when WECS 
generation is present. Thus minimizing WECS genera- 
tion change by WTG model and siting configuration 
selection can either dramatically decrease or possi- 
bly eliminate the need for modifying the unit 
commitment, regulation, and economic dispatch when 
WECS generation is present, 

A brief discussion of spinning reserve, unload- 
able generation, and load following requirements and 
how they are affected by the magnitudes of the change 
and rate of change from an echelon and a farm is now 
presented. The two factors that determine the 
spinning reserve, unloadable generation, and load 
following requirements are the maximum probable rise 

aLj^ and drop ALj^ in thermal load in ten minutes. 
These maximum probable changes are defined as 

< = - 4)T + Op + (4 + Ql^. 0} 

03) 

AL- = max{-(L^^^ - 4)1 + 0^. + 0} 

where 


Q[k Maximum probable drop in system load 

in 10 minutes 

Dp Largest single resource (generation 

or import) subject to failure 

Dq Largest single commitment (export) 

subject to failure 

The spinning reserve $Rj^, unloadable generation 
SSj^, and load following capabilities LF^^ for a unit 
commitment where 

G|^ Required load-following capacity 

f Average minimum generation level of 

load-following units as a fraction of 
maximum capacity 

g. Average operating level of load- 

^ following units above level f(0£g. i 

1 - f) 

r Average ramp rate of load-fol lowing 

units in % of rated capacity per 
minute 

are defined as 

SR, = (1 - f - g,) G, 

LFk = lOr G, 

when there Is no unconnected hydro or pumped storage 
units, interruptible load, and unused but connected 
base loaded generation to contribute to spinning 


reserve and unloadable generation capability. The 
requirements for security on the system are that 

SRk = (i-f-gk)Gk- >AL,; (15) 
SSk = 9kGk±^Lk 06) 
LF, = lOr G, >: max{AL|^, AL”} (17) 
The presence of significant wind generation can 


affect the maximum first contingency loss of resource 
Dp and commitment and the probable rise and 

drop OjJk 1n wind generation in ten minutes. The farm 

penetration constraint [3] limits the maximum change 
in wind generation in an array for passage of a 
thunderstorm to be less than the minimum of the 
maximum first contingency loss of resource or commit- 
ment 


Lk - thermal load at the k hour 


m1n(0p, Dq) 


(18) 
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3. EFFECTS OF THE WTG MODEL 


‘Where is the capacity of each of the wind 

turbines in the farm. If the farm penetration con- 
straint^ Is; vlglited In some region due to favorable 
.economics and limited siting availability due to 
wind, environmental or other factors, then either Dj^ 

or Dq or both must be Increased to the farm capacity 

and Q-^i^ = = 0. The logic for changing Dp, or 

both is that the maximum first contingecny loss of 
resource or conmitment 1$ now the worst case changes 
in wind generation due to passage of the thunder- 
s to rmfs, trail ing edge Inflow and leading edge outflow 
res pectively as Indicated by th^ simulation results 
:r given earlier in this section. and are set to 

zero because the effect of wind generation has been 
already included In adjustment of Dp and D^ and thus 

aL^ and ALj^. 


If the farm penetration constraint is not vio- 
lated, Dp and are not changed and values of and 
Qy^must be determined. These probable or predicted 
changes in wind generation Q^l^) depend on the 

anticipated wind conditions, the wind turbine models 
In the array, the siting of wind turbines, and the 
correlation of wind speeds at the various sites. 
These factors and their effect on the probable change 
In WECS generation in ten minutes will be discussed 
1n the next two sections. 


The analysis In the previous section showed that 
the change and rate of change of wind generation out 
of an echelon and farm are dependent on” the following 
wind turbine model parameters; capacity C^, the slope 


of the change of power produced with wind speed 
and the distance D Into the leading edge 

R Cl 

outflow where the wind speed first reaches rated 
velocity Vp for the WTG model. The capacity and 

distance D are much larger for a MOD-2 WTG than for a 

Cy 

MOD-1 while the slope -n — n — is nearly identical for 


the two WTG models. Thus, the change and rate of 
change of generation for an echelon (10,11) and farm 
(12) of fOD-2 WTGs will be much larger than for an 
identical echelon and farm of MOD-1 WTGSr 


A detailed model and simulation of the MOD-2 and 
the MOD-2 wind farm was developed and is discussed in 
[3], The discussion of the detailed operation of the 
MOD-2 is omitted here. The simulation of the coastal 
farm of MOD-2 wind turbines experiencing the Mitchell 
storm front is presented and compared with similar 
results for the MOD-1 presented in section 2, These 
simulation results on the MOD-2 confirm the results 
of the above analysis of the differences in power 
change and rate of change from the MOD-1 and MOD-2. 


Two methods for determining or setting 

were discussed in [3]. The first method 1s based on 
re lia bil ity analysis that includes the statistics of 
WECS generation change on a particular array for 
anticipated wind conditions and can select in 

addit ion to aL|^, A second method would select Q|Jj^ and 

based on a weighted prediction of change of WECS 

generat ion 'change "on a specffic array and anticipated 
wind conditiorfS, The weighting would depend on the 
operating procedure of the utility. 


!he proper selection of wind turbine models and 
the siting configuration can dramatically reduce 

when the farm penetration constraint is sat- 
1 sf le d; Thus the analysis of how the factors in the 


wind turb i n e nio de 1 ( Cy , 


W’ 1 — siting 
R Cl 


conf Iguration (p^, p^, d) can be selected to reduce 

WECS generation change can provide guidelines for 
generators and model selection. 


The values of Dp and D 

not be reduced by wind turbine model or siting con- 
figuration selection if the farm penetration con- 
straint is violated since the Np ~ ^nd the 

farm penetration level does not necessarily depend on 
the wind turbine mo del and siting configuration but 
ori' the total capacity of wind turbines in an area 
swept out by a thunderstorm. The effective farm 
penetration can be reduced if the farm penetration 
constraint Is violated by the controls discussed in 
section 5, These controls would reduce Dp and D^ and 

thus spinning reserve, unloadable generation and load 
following requirements if the farm penetration con- 
straint mre violated. 


(and thus and iQcan 


The power out of a coastal farm of MOD-2 wind 
turbines experiencing the Mitchell storm front is 
shown in Figure 2, Note that the power increases are 

dp; 

-g| = 1.75 MW/sec. 

which is larger than for the MOD-1 and that the period 
over which power changes on each echelon Is T^ = 80s 

rather than 50s for the MOD-1. Thus the total power 
change, the rate of change, and the period over which 
the change in power on an echelon of MOD-2 WTGs is 
significantly larger than for an identical ehcelon of 
MOD-1 WTGs experiencing the same Mitchell storm front 
as predicted from the analysis. 



Figure 2 The Power Out of a Coastal Farm of 
MOD-2 Wind Turbines Experiencing 
The Mitchell Storm Front 


The power output of the MOD-2 WTGs 1n the first 
echelon fall with the decrease in wind speed after 
passage of the very high winds of the leading edge 
outflow. The power output of a MOD-1 remains con- 
stant during this period because Vp * 18 mph is much 

lower on a MOD-1 causing the blade pitch controls to 
keep power output constant when wind speeds are above 
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The power increase out of the second echelon 
starts at * 220s and is similar to that out of the 

first echelon. The power drop due to shutdown of the 
two echelons is quite rapid. The power decrease out 
of the first echelon is shown at t - 650 seconds 
followed almost immediately by a ramp increase in 
power on the second echelon due to increasing wind 
speeds for passage of the trailing edge inflow over 
this second echelon. 

4. EFFECTS OF SITING CONFIGURATION 
ON WECS GENERATION CHANGE 

The purpose of this section is to discuss how 
following factors affect Jthe probable WECS generation 
change in ten minutes (Qyj^, 

(1) density of WTGs in a farm when N^C^ is 
held constant; 

(2) uniformity of the farm siting configuration 
when N^C^ and values are held constant; 

(3) the wind speed characteristics of typical 
thunderstorms; and 

(4) the correlation of wind speed characteris- 
tics at various sites in an echelon and 
farm. 

This discussion will indicate how each of these 
factors affect probable changes in WECS generation 

for ten-minute intervals and how one 

might attempt to minimize these probable changes and 
thus to minimize the changes in spinning reserve, 
unloadable generation, and load following capability 
required if the farm penetration constraint is satis- 
fied. The siting configuration has no affect on the 
spinning reserve, unloadable generation, and load 
following requirements if the farm penetration con- 
straint (18) is violated since the change in these 
requirements are embodied in changes in Oj^ and to 

the farm capacity (N^C^) and not changes in 
and 

The effects of Increasing and not N^ or 
penetration N^Cy^ is indicated by simulating the WECS 

generation change out of a coastal farm of MOO-1 WTGs 
experiencing a Mitchell storm front when the distance 
between two echelons is decreased from 2 miles to 0.5 
miles. The power change from this modified coastal 
farm, where the total number of WTGs is unchanged but 
the farm density is increased 2.5 times, is shown in 
Figure 3A. The average WECS power change as a 
function of the length of the Interval over v^ich the 
average is computed Is plotted for the WECS genera- 
tion change out of the modified coastal farm siting 
configuration in Figure 3B, Note that the average 
power system response rate capability curve Is also 
included in Figure 3. The Increased farm density has 
not increased the instantaneous rate of change in 
power from an array but held it at the level for 100s 
rather than 50s. The farm and echelon penetrations 
have tq be reduced from 6% and 3>t, respectively, to 4% 
and 2%f respectively, so that the average WECS gen- 
eration change would not exceed power system response 
rate capability. If the density were increased as 
indicated, while holding the number of WTGs (N^) and 

penetration (N^Cy) constant, the peak frequency and 

area control error deviations would be larger and 
would not be reduced to low values as quickly because 


the system response and response rate capability are 
more severely stressed by the same total WECS genera- 
tion change occuring in a shorter interval. 



TIME IIM SECOIVins 



Figure 3 Power Change and Rate of Change 
from the Modified Coastal Farm 
of MOD-1 WTGs 


The effects of reducing both farm and echelon 
siting densities without changing farm penetration is 
shown by simulating the WECS generation changes from 
the midwestern farm and coastal farms of MOD-1 WTGs 
as shown in Figures 4 and 5. The midwestern farm has 
the same 100 WTGs in 10 echelons of 10 WTGs each with 
separation between every WTG equal to 0.7 miles. The 
spreading out of generation, by reducing and 

and reducing to values that more closely approach 
‘’f 

1.4 midwestern 

Pp - 

4 coastal 

eliminates the saturation of area control error, 
reduces frequency and area control error deviations, 
and permits the utility to handle the WECS generation 
easily in a manner similar to load changes. 

Thus, reducing farm and echelon density and 
maintaining echelon density at or near farm density 
levels (uniform low density farm siting patterns) 
could dramatically reduce WECS generation changes in 
ten minutes and the need to increase load following. 
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spinning reserve > and unloadable generation require- 
Mnts/ It should be noted wind speed characteris- 
tics, site availability, legal and environmental con- 
straints ^seriously Timlt the ability to select 
uriifQrm lov^de^ity siting configurations. A coastal 
farm, ' Where ' wi nd speeds “ drop as distance 'from the 
ebast increases, is one example where higher non- 
uniform density configurations are likely. 



Figure 4 Effects of Power from a Coastal Farm 
. _ I ~of) System Frequency and ACE 



figure 5 Effects of Power from a Midwestern 
Farm on System Frequency and ACE 


Tiie siting pattern assumed to this point is that 
all WTGs_ arie in straight lines normal to the 

front although the definition of an echelon included 
all generaTion in strip miles wide and D miles long 

in the iireciion of the motion of the storm front. 
The echeTons were separated by a distance d which was 
greater than the MOD-1 or MOD-2 WTG value of D for the 
coastal and midwestern farm configurations. The 
effects of randomly siting WTGs in the farm main- 
taining V 0.5 mi. separation between WTGs, which 
Avoids turbulence and loss of efficiency, was Inves- 
coastal and midwestern farms. This 
random siting has the effect ^ making the density ^ 

of WTGs in every D x area smaller and much closer 

%9 f ds n s 1 ty ^ev e Ts . Thus, this random siting 

Ih o ws th°e ef f ec t of spreading but the siting within 
fh" echelon and the effects of reducing echelon 
density. 

The coastal farm, with a d = 2.0 mile separation 
between echelons, was randomly sited by restricting 
alT WTGs within the two mile strip but maintaining a 

center band' within this strip to satisfy the 
turbuTehce' avoidance constraint. Note from Figures 1 
and 6 for the original and randomly sited coastal 
farm, respectively, that randomly siting WTGs has 
toothed out the power Increases for passage of the 
thunderstorm leading edge outflow so that the rise 1s 
continuous wfth ho intervals where WECS generation 


change has stopped. Random siting has made the drops 
in generation due to passage of the trailing edge 
Inflow almost as smooth and continuous a$ the in- 
creases. This result Indicates that the large almost 
instantaneous drops in generation, that could occur 
due to simultaneous loss of generation on an entire 
echelon. Is not likely since the siting configuration 
is not likely to be perfectly straight lines normal 
to the motion of the thunderstorm front. The large 
frequency changes, which result due to the Inability 
of frequency regulation to cope with such large 
instantaneous change, 1s also not likely. 

The average response rate over a 70 second 
interval for this random sited coastal farm and the 
original farm are almost identical, as can be seen by 
comparing the power change on the two farms at t = 70s 
in Figures 6 and 1. This indicates random siting has 
virtually no effect on the average rate of change 
that must be coped for by frequency regulation and 
regulation In 60 seconds given by (10) for this 
Mitchell front. However, randomly siting the coastal 
farm would likely have greatly reduced the WECS gen- 
eration changes for the O'Hare 6 wind speed prof 1 le, 
shown In Figure 7B, where adjacent echeTons of the 
coastal farm had simultaneous incre ases. This random 
siting of echelons, which makes echelon densities 
smaller and more equal to farm density, could thus 
have decreased the changes in WECS generation out of 
the coastal farm. This result confirms that 

(a) reducing echelon density and penetration, 
and 

(b) making echelon density uniform by making it 
more nearly equal to farm density 

decreases WECS generation change over ten minutes 

(Qyk’ ^k^ instantaneous rates of change 

which must be handled by frequency regulation and 
regulation controls. 



Figure 6 Random Sited Coastal Farm of 
MOD-2 WTGs with the Mitchell 
Storm Front 

Wind speed time profiles during thunderstorms at 
Chicago's O'Hare International Airport and at 
Ludington, MI were measured and used at Inputs to the 
original coastal farm siting configuration of M0D^2 
_WTGs, replacing the Mitchell storm front wind speed 
profile used exclusively up to this point. The power 
fluctuation out of this coastal farm for these vari- 
ous thunderstorm wind speed versus time profiles are 
given in Figures 7A-7D. 

The actual measurements of wind speeds at these 
sites Indicate there can be several peaks and lulls 
1n a wind profile (Figures 7A-7C ) and that the 
time Interval between peaks can vary between 10 min- 
utes to 40 minutes. The gradual buildup of wind speed 
for an advancing storm and the eventual peaks and 
lulls in the actual storm can be seen In Figure 7D. 



The power from the coastal farm Indicates that 
WECS generation changes can occur simultaneously on 
the two echelons which was not true for the Mitchell 
front. These WECS generation changes sometimes add 
giving short term ( “ 1 minute) WECS generation change 
that is larger than can occur on a single echelon. 
This occurs at t = 300s and 1800s in Figure 7B where 
WECS generation changes reach 150 MW and 170 MW when 
the capacity of the echelon of M0D-2s is 120 MW. The 
repetition or these changes in wind speed and genera- 
tion and the long duration of the thunderstorms (- 1 
hour) were not anticipated based on the Mitchell 
storm front data. 



Figure 7A O’Hare 2 Storm Front on the 
Coastal Farm of MOD-2 WTGs 



Figure 7B O’Hare 6 Storm Front on the 
Coastal Farm of MOD-2 WTGs 



Figure 7C O'Hare 5 Storm Front on the 
Coastal Farm of MOD-2 WTGs 



The effects of these large and cyclic power 
variations from the coastal farm of M0D-2s experi- 
encing the O' Hare 6 wind speed versus time profile is 
shown in Figure 8. The simulation is performed on 
the 4000 MW system with 5% load following capability 
and experiencing i^/min. load change for ten minutes 
in addition to these WECS generation changes on the 
5.55» penetration coastal farm. The power variations 
from the coastal farm with O' Hare 6 wind speed pro- 
file is large and oscillatory. This Is observed in 

large area control error and frequency deviations 
that reach peak to peak 150 MW and .03 hz respectively 
at t = 300s and approximately similar values at 
1800s. The area control error saturated in both 
positive and negative directions within 100s in each 
case. Economic dispatch and economic dispatch/regula- 
tion units take on the load increase over the first 20 
minutes and then respond to the overall cyclic (t = 
800s) power changes in WECS generation, but not the 
faster changes seen on base and hydro units. 



Figure 8 Effect of the O'Hare 6 Storm Induced 
WECS Generation Change on System 
Frequency and ACE 

The large power changes from the coastal farm 
with O'Hare 6 wind speed profile are truly excessive 
for the 4000 MW system not only in_ size but also in 
terms of their repetition at t = 300, 1800 and 2700s 
and the duration of these changes (t = 3600). The 
size of these oscillations is due,^ in part, to the 
occasional overlapping of generation increases or 
decreases on different echelons. These large WECS 
generation changes can repeatedly cause saturation in 
area control error In both directions over a very 
short interval. 

It would appear that reducing farm and echelon 
density without changing farm penetration would re- 
duce the magnitude of these fluctuations. This is 
shown in Figure 9 where the O'Hara 6 wind speed 
profile is inserted into the midwestern farm configu- 
ration. Note that compared to coastal farm changes 
in WECS generation, shown in Figure 9b, the peaks and 
sharp valleys have been eliminated. The result in- 
dicates high echelon and farm density can have a 
major effect on the WECS generation changes over ten 
minutes as well as those over 60 seconds thus in- 
creasing spinning reserve, unloadable generation, and 
load following requirements as well as effecting the 
measures of operating reliability such as the average 
area control error and interval between area control 
error zero crossings. 

It should be noted that the width of peaks and 
valleys, and duration between such peaks and valleys, 
probably are related to the very structure of a 
thunderstorm wind pattern. If statistics were deter- 
mined on the width of peaks and valleys and duration 
between them, rules or principles for siting in 
coastal and midwestern farms could be developed that 
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would nintmize WECS generation change out of an 
array. The analysis of wind speeds and the appro- 
priate principles for siting WTGs is a subject for 
further research. 



Hgure 9 0‘Hare Storm Front on the Midwestern 
0^ MQP-2. WTGS 


Tte and simulation of thunderstorm 
Induce^ wffi S generation changes have assumed that the 
Mjld at every point along a straight line 
ri5rmal^oJ: he dir^ctj on of front motion is identical 
and “perTectly correlated and that the wind speed 
profile j^op agates at so that each WTG observes 

the $ ame w j n d s peed profile. This may not be true in 
tea 1 1 ty an d th us an effort i s made to _ as se s s wh at 
effect th ^ assumption of perfect correlation of wind 
gpeg d in_ an e^el^n and the assumption that the wind 

fi^om echelon to echelon 
Cnchanged has on the size of power variations from a 
farm. The present correlated echelon wind farm model 


works as follows. 



power out of a single WTG Is multiplied by 
In an echelon ^nd then this echelon 

iftjteFBuTpur IS delayed by to get the output 

//I 0 

.of jyi€_ nth echelon. The output of all echelons is 
then siinfjied. This wind farm model assumed all wind 
_^P^^ds In an echelon are perfectly correlated and 
■ Jha t the wind speed proTile propagates from echelon 
t? . A lec ond wind farm model assumes power 

gf WTJ jod^p ood ent. If power out of each 

KT? i s apumg d to be a gaussian process, then if all 
wind^ turbfnes in a ferm see independent identical 
stationary ergodic wind speed processes, a sample 
fyncti^ of the power out of such a farm is 

(Py(t) - (19) 

-1‘ t 

where - Y/pP^(t)dt and P|^(t) is a sample function 

of the Jcwer^^t^W a single WTG for this wind speed 
process. ^ perfectly correJate<f ipqd^ of power out 
o^a f am assumed the wind speed process at every WTG 
are an Tdantical stationary ergodic gaussian 
pr ocesse s are perfectly correlated so that 

( 20 ) 

ihe Output of indopendent inodel, porfoctly 
correlated rrodol, and the perfectly correlated 
l£u|]Sr' Slfjd farm model are shown for wind speed 
profiles measured during thunderstorms at O'Hare 5 
?23 5 ^in fJgures 10 A-C, respectively and 11 A-C. 
Ef?P®i|-lvely. Note that power variations out of the 
perfectly correlated _ model, perfectly correlated 
echelon model , and independent model are generally 
successively smaller. There are exceptions when the 
perfectly correlated echelon model has larger power 
variations than the perfectly correlated model, which 
occurs when power increases on different echelons 
simultaneously. 


The assumptions concerning the correlation of 
wind speeds at various sites thus not only affect the 
magnitude of the variations over a one minute inter- 
val but also the variations over a ten minute inter- 
val. ITius, the correlation of wind speeds can siq- 
nificantly affect the statistics of the changes in 

WECS generation out of an array over a ten minute 
interval and thus the selection of Q"|^) in the 

load following, spinning reserve, and unloadable gen- 
eration requirements. 



Figure 10 Independent, Correlated Farm, and 
Correlated Echelon Wind Speeds on 
WECS Generation from the Coastal 
Farm with O'Hare 5 Wind Speed Profile 



Figure 11 Independent, Correlated Farm, and 
Correlated Echelon Wind Speeds on 
WECS Generation from the Coastal 
Farm with O'Hare 6 Wind Speed Data 


5. EFFECT OF WIND ARRAY CONTROLS 


me previous two sections discussed factors, 

which depend on the WTG model (C^, 0. .) and 

\ W Vr - Vci 

siting configuration (p^, p^. d). that affect prob- 
able WECS generation change The selec- 

turbine models and siting configuration 
3^sed on economics, wind conditions, site 
availabiltiy and other factors which do not permit 
.WTG selection and siting 
configuration combination. Moreover as WECS genera- 

configuration will reduce the probable changes (^^. 

Q^l^) in WECS generation sufficiently to eliminate the 

reserve, unloadable gen- 
eration and load following requirements and the 
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appropriate modification of regulation and economic 
dispatch controls as described In [3]. Finally, WTG 
model and siting configuration selection only have 
effect on reducing the need for modifying unit 
commitment, regulation, and economic dispatch if the 
farm penetration constraint is satisfied. 


The direct controls of WECS generation change, 
discussed in this section, can: 

(1) reduce the effective farm penetration when 
thunderstorms are present and thereby make 
an array that would otherwise violate the 
farm penetration constraint effectively 
satisfy the constraint. This satisfaction 
of the farm penetration constraint thus 
makes the Increase in spinning reserve, 
unloadable generation, and load following 
requirement depend on the probable ^hange 
in WECS generation in ten minutes and 

and not on the nwdif ication of Dp and 
Dc to Np C^; 

(2) significantly reduce the probable W|CS gen- 
eration changes in ten minutes (Qyj^» Qyj^) 

assuming that WTG model and siting configu- 
ration have been appropriately selected and 
that the effective penetration of the array 
with these controls present satisfies the 
farm penetration constraint. 

Two direct controls limit the WECS generation 
change out of an array during any ten minute interval 
by coordinated control of blade pitch angles of all 
WTGs in an array. These controls would also reduce 
the apparent farm penetration during thunderstorms so 
that it meets the farm penetration constraint. One 
of these controls would clip WECS generation change 
in any ten minute interval and the other causes 
partial shutdown of each echelon in anticipation of 
the storm so that the effective farm penetration 
satisfies the farm penetration constraint and the 
WECS generation changes are capable of being handled 
by the unit commitment, regulation, and economic dis- 
patch controls that are set when WECS generation is 
not present. 
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Figure 12 Effect of Partial Shutdown on 
ACE and Frequency Deviation 


6. CONCLUSIONS 

The paper discusses methods of reducing the wind 
generation changes from an array for passage of a 
thunderstorm by wind turbine model selection and site 
configuration selection. Coordinated blade pitch 
controls are also discussed and can be used to reduce 
the effective farm penetration level so that the farm 
penetration constraint is not violated. Satisfaction 
of this constraint implies that the added spinning 
reserve unloadable generayon and load following re- 
quirement on unit commitment and the added response 
capability of AGC controls depends on WTG model se- 
lection and site configuration selection and not on 
the capacity of the farm. These coordinated blade 
pitch controls on each WTG in the array could also 
reduce the wind generation change out of an array 
much as wind turbine model and site configuration 
selection. The wind generation change after appro- 
priate wind turbine model selection, site configura- 
tion selection, and coordinated blade pitch controls 
must be responded to by the units under AGC control. 
Limitation on site availability and wind turbine 
model selection and economic incentives for higher 
density siting may contribute to a rather significant 
change and rate of change in WECS generation espe- 
cially during severe weather conditions. The adjust- 
ment of unit commitment to allow sufficient spinning 
reserve, unloadable generation, and load following 
capability and the adjustment of AGC controls to 
exploit the response capability available from unit 
commitment are discussed in [3J. 


The partial shutdown of each echelon reduces 
both the ramp WECS generation increases and the 
sudden WECS generation drops on each echelon by 50%. 
The area control error, frequency, and tie line power 
deviations for each of these changes is thus 
decreased by approximately 505! also. 

Frequency regulation is seen to be more capable 
of quickly reducing frequency deviations after the 
sudden drops of WECS generation on each echelon be- 
cause they are smaller. Finally the saturation of 
area control error after the drop in WECS generation 
on the second echelon is reduced from 500s to 100s 
indicating the effectiveness of a partial shutdown in 
anticipation of the arrival of a thunderstorm front. 

This partial shutdown would require wind speed 
monitors .to detect the approach of a thunderstorm 
from any direction. Both clipping and partial shut- 
down would result in lost energy and that may some- 
what reduce the economic attractiveness of the wind 
turbine arrays. 
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